Spin-polarized tunneling current through a thin film of a topological insulator in a 

parallel magnetic field 
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We calculate the tunneling conductance between the surface states on the opposite sides of an 
ultra-thin film of a topological insulator in a parallel magnetic field. The parallel magnetic field 
produces a relative shift of the in-plane momenta of the two surfaces states. An overlap between 
the shifted Fermi circles and spinor wave functions result in unusual dependence of the tunneling 
conductance on the magnetic field. Because spin orientation of the electronic surface states in 
topological insulators is locked to momentum, spin polarization of the tunneling current can be 
controlled by the magnetic field. 
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I. INTRODUCTION. 
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Iii the past few years, there has been a rapid 
progress in the field of topological insulators (TIs) 
The three-dimensional (3D) topological insulators have 
the helically-spin-polarized surface states with the two- 
dimensional (2D) Dirac dispersion, which are observed 
experimentally [2j. TIs not only offer an exciting play- 
ground for fundamental physics, but also have a variety 
of potential applications ranging from spintronics [3l4l2j 
to quantum computing [1 31 ] - However, the experiments 
show that the Fermi level is often lifted to the conduc- 
tion band, thus making the topological surface states less 
relevant for the properties of real materials [Tij . In part, 
this motivated the study of the thin films of TIs, where 
the Fermi level, bulk gap and the hybridization between 
the opposite surface states can be tuned [15M171 ] . A num- 
ber of intriguing effects for the TI films were predicted 
in the case where the top and bottom surfaces interact 
[l8l ]. In the ultra-thin limit, a gap opens in the surface 
states electronic spectrum due to hybridization between 
the opposite surfaces [lj| [l9l - [2lT | . 

While most papers study the in-plane transport prop- 
erties of the surface states, here we study the tunneling 
conductance G per unit area across the ultra-thin TI film 
in the presence of a parallel magnetic field B y (see Fig.[l|. 
The in-plane magnetic field results in a relative shift of 
the Dirac cones of the surface states in momentum space, 
as discussed in related works on graphene bilayer [22.123$. 
However, unlike the spectrum of graphene, the spectrum 
of the TI thin film remains gapped until a critical value 
of the magnetic field is achieved [24J . Motivated by the 
original measurements of G vs B y for a GaAs bilayer 
by Eisenstein et al. [25j], here we calculate the tunnel- 
ing conductance for a TI film as a function of the par- 
allel magnetic field. In contrast to other systems, the 
spin structure of the electronic spectrum for the surface 
states of TI results in unusual dependence of the tunnel- 
ing conductance on B y . In a recent experiment (26[, the 
authors observed a strong B y dependence of the out-of- 
plane current in a film of Bi2Te3. The predicted theoret- 
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FIG. 1: (Color online) Thin film of a topological insulator of 
thickness d in a parallel magnetic field B y . The surface states 
tpi and -02 shown in red and blue overlap and couple when 
the thickness d is comparable with the decay length £ of the 
surface state. 



ical dependence G vs B y is in qualitative agreement with 
the experimental curve. In addition, we predict that the 
spin polarization of the tunneling current is proportional 
to the magnetic field. The 100% spin polarization can be 
achieved for an experimentally accessible strength of the 
magnetic field. 



II. THE MODEL. 

Let us consider a film of TI, which has two states ip\ 
and tp2 localized at the opposite surfaces of the film, as 
shown in red and blue in Fig. [1] While the wave func- 
tions are localized in the z direction, electrons are free to 
move parallel to the surface, so the in-plane momentum 
p = (p x ,p y ) is a good quantum number. For clarity of 
consideration, we model the surface states by the simple 
Rashba Hamiltonians 

h i>ui>2 = ±h (p), (1) 
h(p) = v z(& x p) = v(& x p y - 6- y p x ). (2) 
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Different signs ± correspond to the surface states ipi and 
ip2 and describe the unit Rashba vectors ±z, which are 
collinear with the normals to the corresponding surfaces. 
The variable v has dimensions of velocity, and a a = (1, <x) 
denotes a full set of the operators acting in the spin space, 
where 1 is a 2 x 2 unit matrix and <r = (a x , a y , a z ) are 
the Pauli matrices. Hamiltonian ^ has the following 
eigenstates and eigenenergies 



P,±) = 



( 
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-i arctan(p x /p H ) 



Eo(p) = ±v\p\ 



where the spectrum has a well-known form of the 2D 
Dirac cone. The spinors in Eq. ([3]) describe the helical 
spin-momentum locking with the spin polarization per- 
pendicular to the momentum 
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so the helicity ((3]) is opposite for the | p, +) and | p, — ) 
bands. Because of the minus sign for in Eq. ([TJ, 
the bands | p, ±) are energetically-inverted, and, so, the 
helicities for the opposite surface states ipi and -02 are 
opposite for a fixed energy. The wave functions of the 
surface states tpi and ip2 decay into the bulk and have a 
finite decay length £, as illustrated in Fig. [TJ So, when 
the thickness of the film d becomes comparable with the 
decay length d ~ £, there is a finite coupling between the 
surface states w, which we assume to be proportional to 
the unit matrix in the spin space. So, the TI film can be 
modeled by a 4 x 4 Hamiltonian H(p), which acts on the 
4-component wave function \l/ 



H(p) = ( H{P) 



w 

-Hp) 



f = 



1p2 



(5) 



This Hamiltonian was also used in Ref. [24| for descrip- 
tion of a thin film. Even though the bulk bands usually 
contribute to the electronic properties of the real ma- 
terials [IH, we ignore a possible bulk effect in order to 
highlight the contribution of the surface states. Hamilto- 
nian ([5]) has the spectrum E(p) = ±y v 2 p 2 + w 2 , where 
the energy gap is determined by the tunneling element 
w. Experimentally, gap varies from 0.25 eV for the ultra- 
thin 2 nm film to 0.05 eV for the 5 nm film of Bi 2 Se 3 [U . 



III. SPECTRUM IN A PARALLEL FIELD. 



Let us now discuss the spectrum of the TI thin film 
in a parallel magnetic field B = B y y. In the Landau 



gauge A — zByX, the Peierls substitution p 
transforms the Hamiltonian ([5]) into 



p — eA 



Hip) 



h(p-q/2) 



-hip + q/2) 



, q 



eBydx. 



The Zeeman term can be included in the definition of q 
[24l [27| . Notice that the parallel magnetic field leads to 



E Pj 




FIG. 2: (Color online) Electronic spectrum Q of the thin 
film consists of the two Dirac cones spin-polarized in the op- 
posite directions as shown with circulating arrows. The Dirac 
cones are shifted by q x due to the parallel magnetic field B y . 
Tunneling between the Dirac cones is dominated by the elec- 
trons with momenta pi and P2 where the shifted Fermi circles 
\p ± q/2\ — pf intersect. 



the opposite momentum shifts Ap x = ±q/2 on the op- 
posite surfaces located at z = ±d/2. Indeed, consider an 
electron traveling between the surfaces with the classical 
velocity v z in the z direction. Then, the Lorentz force 
Fl = ev z B y x results in the in-plane momentum change 
q = J dtFL = eBydx between the opposite surfaces. So, 
in the limit w — > 0, the spectrum of the Hamiltonian ((5]) 
is given by the two shifted Dirac cones 



E{p) 



±v 



p± 



(7) 



as shown in Fig. [2j A similar shift of the Dirac cones 
was discussed for graphene multilayer due to a parallel 
magnetic field [22|, [23| and for a twisted graphene bilayer 
[28j . In general, when w is not small, the spectrum of 
Hamiltonian ([6]) was calculated in Ref. [24|. However, we 
focus on the former case when w is small and include w 
as a perturbation. 



IV. TUNNELING CURRENT IN A PARALLEL 
FIELD. 



When a small potential difference V is applied between 
the opposite surfaces of the film, the out-of-plane tunnel- 
ing current per unit area 

3a = (30,3) = iw (e ieV Vi<^V>2 - e- ieVt 4a a ^x) (8) 

flows between the opposite surfaces, where ipi and ^2 are 
the annihilation operators on the corresponding surface 
states. Here, we put h = 1 for simplicity and reestablish 
the correct units in the end of calculation. So, the tun- 
neling current is written in a convenient four-component 
form, such that ejo and -|j are the charge and spin cur- 
rents. Using the linear-response theory for w as a per- 
turbation [29(, the tunneling conductances per unit area 
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FIG. 3: (Color online) (a) The tunneling conductances for 
spin-unpolarized Gq 1 ', Eq. (|12[1 . and spin-polarized Gq 2 ', 
Eq. (|13[) . Fermi circles are shown by the dashed and solid 
black lines, respectively. The latter curve G corresponds 
to the negative magnetoresistance. The red curve is the spin 
polarization (|15p of the tunneling current, and the right ver- 
tical axis shows its value, where 1 corresponds to 100% spin 
polarization, (b) Schematic views of the shifted Fermi circles 
\p ± q/2| = pf for the corresponding values of q. The tun- 
neling is allowed for the momenta pi and p2 where the Fermi 
circles intersect, also shown in Fig. The spin polarizations 
for these momenta are shown by the blue and red arrows cor- 
responding to the ipi and ip2 surfaces. The vector sum of the 
arrows of the same color defines the net spin polarization of 
the tunneling current, which grows with the increase of q as 
illustrated in the bottom of panel (b). 



G a = limy_>o (ja) /V can be expressed through the spec- 
tral functions Ai^(E,p) evaluated at the Fermi energy 
Ep as 



J d 2 P Tr \& a Ai (e f , P + I) A 2 ( 



(9) 

Here, G a = (Go, G) has four components, where eGo and 
■|G correspond to the tunneling conductance for charge 
and spin currents. The spectral functions A±,2 are given 
in the momentum space 



■■ :^. ! p:|)= 5{v\ P t\\-E f ) 

x | P t|,±)(pt|,± 



(10) 



where the opposite signs correspond to the opposite sur- 
faces f/'i a n d ip2i and the spinors | p, ±) are defined in 
Eq. ©. Using Eqs. © and (HU]), let us first calculate the 
tunneling conductance Go for the charge current 



Go — 



2tt r 2 



Jd 2 P S> 

f(p) = 



{p- 



Q 
2 • 



Pf) $ 
\P + 



P+ f | - Pf) f(p) 

v |2 



where p F — Ep/v is the Fermi momentum. Accord- 
ing to Eq. (fnj) . the tunneling current is carried by the 
electrons that have the in-plane momenta p\ and p 2 cor- 
responding to the intersection points of the shifted Fermi 
circles |p ± q/2\ = pp, as shown by the thick dots in 
Fig. [2] In addition, there is a contribution f(p) due to 
the non-trivial scalar product between the spinors ([3]) 
corresponding to the Fermi circles. 

As a warm up, let us first consider the case where is no 
spin- momentum locking, so f(p) = 1 and Eq. (|llj) gives 



2w epp 



Trh 3 v 2 qy /, App — q 2 



(12) 



The tunneling conductance Gq diverges at q = and 
q = 2pF, as shown by the dashed line in Fig. EJa). For 
a small magnetic field q 2pp, the tunneling conduc- 
tance is large because of the large overlap between the 
two Fermi circles shown in the left part of Fig. [3][b) . For 
an intermediate magnetic field q ~ pp, the Fermi circles 
intersect only at two points p\ and P2, so the tunnel- 
ing conductance decreases. For q ~ 2pp, the two points 
Pi and P2 come together, so the Fermi circles overlap 
at the locally flat regions as shown in the right part of 
Fig. IHb), and the tunneling conductance is large again. 
Once q > 2pp, the two Fermi circles separate, and the 
tunneling conductance drops to zero. The experimental 
curve of the tunneling conductance for a GaAs bilayer 
[25l ] shows behavior similar to Eq. (fT2"j) . but the diver- 
gences at q — and q = 2pp are smeared out due to a 
finite scattering time r. 

Now, let us consider the form- factor f(p) in Eq. (TTTj) 
arising from the spinor eigenstates ([3]) of the Rashba 
Hamiltonian ([2|). As discussed above, the surface states 
have opposite spin polarizations. In Figs. [5] and [31 the 
polarizations corresponding to the different Fermi circles 
are shown by the blue and red arrows intersecting at the 
angle (f>, the value of which follows from simple trigonom- 
etry: cos0 = q/2pp. The form-factor f(p) in Eq. (TTTj) 
takes into account an overlap between the spinors on the 
different Fermi circles and is equal to /(pi.2) = cos 2 0. 
So, the tunneling conductance (jTTJ) is multiplied by a 
factor of q 2 /Ap 2 F relative to Eq. (TT! 



G ( 2) (q) 



w 2 eq 



2Trh 3 v 2 ^/Ipp 



(13) 



(11) 



The tunneling conductance given by Eq. (|13j) is plotted 

by the solid line in Fig.[3ja). In contrast to Gq , the tun- 

(2) 

neling conductance G is suppressed as q — > 0, because 
spin polarizations of the Fermi circles are opposite at the 
points pi and P2, as shown in Fig E](b) on the left. This 
is a signature of the opposite spin helicity of the Fermi 
circles and thus may be used as a way to detect it. For 
q near 2pF, however, the spin configurations are almost 
aligned, as shown in Fig. EJb) , so, Gq 1 '' and G Q 2 ^ become 
equal. 
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The non-trivial spin structure of the Fermi circles not 
only modifies the tunneling current, but also makes it 
spin-polarized. Let us define spin polarization as the ra- 
tio of the spin current over the charge current 



where the second equation follows from the definition 
of the conductance. Using the property that the eigen- 
states ([3]) have the well-defined spin orientation we 
obtain from Eq. © 

S(p! - g/2) + S(p 2 - g/2) ed 
S = = cos<^ = — B. 

(15) 

So, the spin polarization of the tunneling current is de- 
termined by the vector sum of the spin directions at the 
momenta p\ and P2- Indeed, in the process of tunnel- 
ing, the electrons with the in-plane momenta p\ and P2 
move from the "blue" to "red" Fermi circle, carrying the 
total spin, which is a vector sum of the spins at the mo- 
menta p\ and P2 , in agreement with Eq. (|15l) . The helical 
spin configuration Q is such that the net spin polariza- 
tion S is parallel and proportional to B, as illustrated in 
Fig-EIb). At q = 0, the spins at the points pi and p 2 are 
opposite, and the spin polarization of the tunneling cur- 
rent vanishes. At q = 2pp , the spins at the points pi and 
Pi are collinear and the tunneling current is fully spin- 
polarized. Also, notice that the electron tunneling at p\ 
changes the spin polarization from the "blue" to "red" 
arrow, thus creating a torque S(p\ + g/2) — S(p\ — g/2) 
in x direction. However, the torque is opposite at P2, so 
the net torque is zero. 



which is corroborated by an increase of resistance R zz 
in the experiment. At the same time, the magnitude of 
negative magnetoresistance increases dramatically with 
doping, as shown in Fig. 4(a,b) of Ref. (26|. This obser- 
vation is consistent with Eq. (|13|) . where the conductance 

(2) 

G increases when pp decreases (for q < 2pp). How- 
ever, our idealized model may be not fully applicable 
to the experiment (26j . where polycrystalline films were 
studied. Grain boundaries and defects may host topolog- 
ical states [U HH , which can contribute considerably to 
the tunneling current. 

Equation (|13|) was obtained for an idealized situation 
where the temperature T and the inverse scattering time 
1/t are much lower than the Fermi energy Ep. For re- 
alistic TI materials, the scattering rate 1/t is finite due 
to impurities or other mechanisms. When 1/t, T, and 
fluctuations of chemical potential [30j become compara- 
ble with Ep, the predicted effects would be smeared out. 
However, the majority of TIs are not strongly-correlated 
systems, so the effects of interactions between electrons 
should not alter the predicted effects considerably. 

While the isotropic Dirac cone approximation is valid 
for small energies in Eq. (|3|), it is not exact for higher 
energies, where energy spectrum has hexagonal warp- 
ing anisotropy, and spin polarization has an out-of-plane 
component. In this case, the analytic formula ([13")) is 
not applicable, but can be easily generalized. The ap- 
propriate spectrum E(p) and spinors | p, ±) should be 
substituted into Eq. (fTO)) , and the tunneling conductance 
obtained from Eq. ©. In general, the tunneling con- 
ductance may depend on the in-plane orientation of the 
magnetic field B due to anisotropy of the surface state 
dispersion. This effect can be utilized to obtain informa- 
tion about hexagonal warping experimentally. 



V. EXPERIMENTAL RELEVANCE. 



Let us estimate the critical magnetic field B y where 
the Fermi circles almost detach, i.e., where q = 2pp. We 
take the realistic value v = 5 x 10 5 m/s, assume d = 5 nm 
(5 quintuple layers of Bi2Se3) and carrier concentration 
corresponding to Ep = 30 mcV. Using Eq. ([6]) for q 7 
we estimate the critical value of the field as B y = 20 T. 
This value is experimentally accessible and can be further 
reduced by either increasing the thickness d or decreasing 
the Fermi momentum pp. 

An intriguing strong negative magnetoresistance effect 
was reported for the Sn-doped films of Bi2Te3 in Ref. [26| . 
A weak in-plane magnetic field less than 1 T causes a 
large drop of the out-of-plane resistance R zz as shown 
in Fig. 4(c) of Ref. [26]. Since R zz oc 1/Go, a decrease 
in resistance R zz corresponds to an increase in conduc- 

(2) 

tance Go, which qualitatively agrees with plot of G 
vs q shown in Fig. [3][a) of our paper. In addition, the 
effect of Sn doping was studied in Ref. (2(| . The Sn dop- 
ing leads to a decrease of carrier concentration and pp, 



VI. CONCLUSIONS. 

We have theoretically studied tunneling conductance 
between the opposite surface states in a thin film of TI 
when a parallel magnetic field is applied. The helical spin 
polarization and the overlap between the Fermi circles re- 
sult in the unusual spin polarization of the tunneling cur- 
rent. Our theoretical results are qualitatively consistent 
with the experiment [2^. However, further experimen- 
tal verifications are needed. First, magnetoresistance in 
Ref. [26[ was measured for a relatively weak magnetic 
field B < 1 T. We predict non-monotonous behavior of 
magnetoresistance for a large magnetic field: the conduc- 
tance (resistance) could sharply decrease (increase) when 
the magnetic field exceed the condition B > 2pp/ed. 
Second, a measurement of the spin polarization of the 
tunneling current is desirable. This effect may be impor- 
tant for spintronic applications and may pave the way to 
observation of spin-polarized currents in TIs. 
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